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1  Main  Theme 


In  many  military  and  biological  applications  of  signal  processing,  one  is  faced  with  the  situation  that 
the  signals  recorded  by  the  sensors  are  largely  influenced  by  how  the  system  of  interest  is  probed  or 
excited.  Therefore,  besides  design  of  optimal  estimators,  it  is  possible  to  acquire  information  about 
the  system  more  efficiently  by  design  of  input  signals  and  waveforms  that  probe  the  system  of  interest. 
These  input  signals  should  be  so  designed  that  they  attempt  to  reduce  the  uncertainty  of  the  unknown 
parameters  the  most.  In  this  proposed  work,  we  aim  to  develop  a  mathematical  theory  which  will 
lead  to  optimal  probing  procedures  with  applications  ranging  from  military  applications  including 
radar  technology  to  biological  applications  like  protein  NMR  spectroscopy.  Since  the  inception  of 
the  grant,  our  main  focus  has  been  to  study  design  of  input  probe  signals  (pulse  sequences)  in  NMR 
spectroscopy  that  are  robust  to  inhomogeneity  and  dispersion  in  the  parameters  of  the  spin  systems 
that  we  are  probing. 


2  Applications 

2.1  Robust  NMR  Spectroscopy 

In  NMR  spectroscopy,  probe  signals  (radio  frequency  pulses)  are  used  to  manipulate  the  dynamics 
of  the  nuclear  spins.  The  goal  is  to  disturb  the  system  from  its  equilibrium  position  and  by  observ¬ 
ing  the  return  dynamics  of  the  system  one  can  infer  something  about  the  system  parameters.  A 
major  challenge  in  these  problems  is  to  design  a  probe  signal  that  produces  a  uniform  excitation 
of  an  inhomogeneous  ensemble  of  systems  with  large  dispersion  in  their  parameters.  For  example, 
in  magnetic  resonance  experiments,  the  spins  of  an  ensemble  may  have  large  dispersion  in  their 
natural  frequencies  (Larmor  dispersion),  strength  of  applied  rf-field  (rf-inhomogeneity)  and  the  re¬ 
laxation  rates  of  the  spins.  In  solid  state  NMR  spectroscopy  of  powders,  the  random  distribution  of 
orientations  of  inter-nuclear  vectors  of  coupled  spins  within  an  ensemble  leads  to  a  distribution  of 
coupling  strengths.  A  canonical  problem  in  control  of  quantum  ensembles  and  quantum  information 
processing  is  to  develop  external  excitations  that  can  simultaneously  steer  the  ensemble  of  systems 
with  variation  in  their  internal  parameters  from  an  initial  state  to  a  desired  final  state.  From  the 
standpoint  of  mathematical  control  theory,  the  challenge  is  to  simultaneously  steer  a  continuum  of 
systems  between  points  of  interest  with  the  same  control  signal.  In  many  cases  of  practical  interest, 
one  wants  to  find  a  control  field  that  prepares  the  final  state  as  some  desired  function  of  the  pa¬ 
rameter.  For  example,  slice  selective  excitation  and  inversion  pulses  in  magnetic  resonance  imaging 
The  problem  of  designing  excitations  that  can  compensate  for  dispersion  in  the  dynamics  is  a  very 
important  problem  in  NMR  spectroscopy  and  is  of  fundamental  importance  in  developing  a  theory 
of  optimal  probe  signal  design. 

In  last  couple  of  years  we  have  had  a  major  conceptual  breakthrough  in  understanding  of  what 
aspects  of  system  dynamics  makes  it  possible  to  design  probe  signals  signal  that  can  compensate  for 
dispersion  in  the  system  dynamics.  We  have  used  these  ideas  for  design  of  pulse  sequences  in  both 
liquid  and  solid  state  NMR  spectroscopy  which  give  far  superior  performance  than  state  of  the  art 
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methods  [6-14].  This  has  resulted  in  over  two  dozen  publications  in  refereed  journals  over  the  period 
of  last  2  years.  We  sketch  the  main  ideas  of  our  work  in  the  following  section.  Then  we  show  how 
we  have  used  these  ideas  in  developing  methods  for  NMR  and  MRI  in  Inhomogeneous 
fields  which  is  the  major  breakthrough  in  our  research  in  the  past  year. 


3  Main  Ideas:  Synthesis  of  higher  powers  of  dispersion  pa¬ 
rameters  by  Lie  Bracketing 


To  fix  ideas,  consider  the  following  equations  that  model  evolution  of  spin  ^  in  a  magnetic  field.  Let 
X  denote  the  unit  vector  pointing  in  the  direction  of  magnetization  of  the  spin,  and  v,u  represent 
strength  of  applied  magnetic  field  along  x  and  y  axis.  In  practice  the  magnetic  field  applied  to  the 
system  is  not  uniform  across  the  sample  and  therefore  different  spins  see  a  different  magnetic  field 
(rf-inhomogeneity)  modeled  by  a  scale  parameter  e. 

X  =  e(u(t)ny  -f  v(t) Flx)X 

where 
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are  the  generators  of  rotation  around  x ,  y  and  2  axis,  respectively.  Consider  now  the  problem 
of  designing  controls  u(t)  and  v(t)  that  simultaneously  steer  an  ensemble  of  such  systems  with 
dispersion  in  their  natural  frequency  and  strength  of  rf-field  from  an  initial  state  X(0)  =  (0,0, 1)  to 
a  final  state  Xp  =  (1,0,0). 

Observe  for  small  dt,  the  evolution 

U\(\fdt)  =  exp(-efly\/di)  exp (-eQx\/di)  exp(efly  \fdt)  exp(eQx\/di) 

to  leading  order  in  dt  is  given  by  I  -I-  (dt)[eQy ,  eflx],  i.e.,  we  can  synthesize  the  generator  [eilx,  efly]  = 
e2fl2,  by  back  and  forth  maneuver  in  the  directly  accessible  directions  Qx  and  Qy. 

Similarly,  the  leading  order  term  in  the  evolution 

U2  =  U\{  —  \fdt)  exp(— eflydt)Ui(\/di)  exp(eQydt). 

is  [eQy,  [eflx,£fh/]]  =  e3flx.  Therefore  by  successive  Lie  brackets,  we  can  synthesize  terms  of  the  type 
e2fc"l"10I.  Now  using  {eflx,e3flx, . . .  ,€2n+1flx}  as  generators,  we  can  produce  an  evolution 

exp{^c/t€2fc+1nx}, 
fc= 0 

where  n  and  the  coefficients  c*  can  be  chosen  so  that 

f^cke2k+1^e  (1) 

k= 0 

for  all  e  €  [1  —  S,  1  +6].  Hence  we  can  generate  an  evolution  exp(0flx)  for  all  e  to  any  desired  accuracy. 
Therefore,  we  achieve  robustness  with  dispersion  to  e  by  generating  suitable  Lie  brackets.  Similar 
arguments  show  that  we  can  generate  any  evolution  exp (/3fly)  and  as  a  result  any  three  dimensional 
rotation  in  a  robust  way.  It  is  also  now  easy  to  see  that  we  can  synthesize  rotation  0  with  a  desired 
functional  dependency  on  the  parameter  e.  Parametrize  a  rotation  in  ©  €  50(3)  by  the  Euler  angles 
(a,/?, 7)  such  that  0  =  exp(aftx)  exp(/3Qy)  exp(7flx).  Given  continuous  function  (a(£), /?(e),7(e)) 
of  £,  we  can  find  polynomials  that  approximate  a(e),  /3(e)  and  7(e)  arbitrarily  well  and  use  these  to 
generate  a  desired  rotation  0(e)  as  a  function  of  e.  Hence  there  exists  a  control  field  that  maps  a 
smooth  initial  distribution  X0(e)  to  a  target  distribution  Xp(e). 
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4  Mobile  Magnetic  Resonance  Imaging  and  Sensing 


Modern  magnetic  resonance  imaging  and  spectroscopy  systems  typically  demand  the  placement  of 
the  sample  or  subject  inside  the  bore  of  a  large  superconducting  magnet  producing  a  homogeneous 
magnetic  field.  This  requirement  makes  high-field  magnets  bulky  and  expensive.  For  many  applica¬ 
tions,  it  would  be  useful  if  a  mobile  magnet  could  be  scanned  over  an  otherwise  inaccessible  object 
in  order  to  acquire  magnetic  resonance  information.  The  advantage  of  such  ex-si tu  analysis  is  that 
limitations  in  sample  size  and  portability  no  longer  prevail.  Mobile  NMR/MRI  systems  essentially 
about  bringing  the  magnet  to  the  object  rather  than  bringing  the  object  to  the  magnet.  Such  a 
system  will  serve  multiple  roles  in  a  combat  situation  including  diagnosing  head  injury  and  trauma, 
detection  of  explosives  and  scanning  metabolites  in  the  living  tissue. 

One  principal  problem  with  portable  or  one-sided  magnets  is  that  the  magnetic  fields  are  neces¬ 
sarily  spatially  inhomogeneous.  As  a  consequence,  the  NMR  spectra  become  broadened  to  the  extent 
that  resolution  is  lost  and  the  chemical  shift  information  is  hidden.  We  are  developing  new  methods 
of  controlling  spin  dynamics  that  allows  the  observation  of  resolved  NMR  spectra  and  MRI  pictures 
of  an  object  or  subject  in  the  inhomogeneous  field.  Our  work  focuses  on  pulse  and  field  schemes  that 
enable  the  retrieval  of  spectroscopic  information  even  when  the  static  field  (BO)  and/or  the  rf  field 
(Bl)  are  strongly  inhomogeneous.  These  methods  rely  on  the  fact  that  the  response  of  the  spins  to 
driving  excitation  is  nonlinear  which  helps  to  develop  external  excitations  that  compensate  for  their 
own  inhomogeneities. 

In  past  year  we  have  developed  control  algorithms  based  on  Fourier  synthesis  for  designing  phase 
compensating  radio-frequency  (rf)  pulse  sequences  for  high-resolution  nuclear  magnetic  resonance 
(NMR)  spectroscopy  in  an  inhomogeneous  B0  field.  We  show  that  using  radio  frequency  pulses  and 
time  varying  linear  gradients  in  three  dimensions,  it  is  possible  to  impart  the  transverse  magnetiza¬ 
tion  of  spins  a  phase,  which  is  a  desired  function  of  the  spatial  (x,y,z)  location.  Such  a  sequence 
can  be  used  to  precompensate  the  phase  that  will  be  acquired  by  spins  at  different  spatial  locations 
due  to  inhomogeneous  magnetic  fields.  With  this  precompensation,  the  chemical  shift  information 
of  the  spins  can  be  reliably  extracted  and  high  resolution  NMR  spectrum  can  be  obtained.  This 
makes  NMR  in  inhomogeneous  fields  possible  [2]. 

In  next  3  years,  we  will  develop  a  light  weight  mobile  magnetic  resonance  imaging  and  spec¬ 
troscopy  system  suitable  for  battle  field  operation.  Such  a  imaging  system  with  help  diagnose 
mental  trauma,  injury  and  fracture  in  a  battle  field  scenario.  Magnetic  resonance  sensing  will  be 
used  for  detection  of  explosives  and  chemicals  and  analysis  of  metabolites  in  the  living  tissue.  The 
enabling  technology  for  developing  such  a  system  is  novel  control  techniques  for  manipulation  of 
spin  dynamics  in  inhomogeneous  fields  in  such  portable  imaging  systems.  These  control  techniques 
exploit  in  a  fundamental  way  noncommutativity  and  nonlinearity  in  the  dynamics  of  spins  as  a 
function  of  external  excitations  [1,  2]. 
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